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ABSTRACT 

We calculate many different nova light curves for a variety of white dwarf masses and chemical compositions, 
with the assumption that free-free emission from optically thin ejecta dominates the continuum flux. We show 
that all these light curves are homologous and a universal law can be derived by introducing a "time scaling 
factor." The template light curve for the universal law has a slope of the flux, F oc f" 1 75 , in the middle part 
(from ~ 2 to ~ 6 mag below the optical maximum) but it declines more steeply, F oc f~ 3,5 , in the later part 
(from ~ 6 to ~ 10 mag), where t is the time from the outburst in units of days. This break on the light curve 
is due to a quick decrease in the wind mass loss rate. The nova evolutions are approximately scaled by the 
time of break. Once the time of break is observationally determined, we can derive the period of a UV burst 
phase, the duration of optically thick wind phase, and the turnoff date of hydrogen shell-burning. An empirical 
observational formula, t-$ = (1.68 ±0.08) ?2 + (1.9 ± 1.5) days, is derived from the relation of F oc t , where ti 
and ?3 are the times in days during which a nova decays by 2 and 3 mag from the optical maximum, respectively. 
We have applied our template light curve model to the three well-observed novae, VI 500 Cyg, VI 668 Cyg, 
and V1974 Cyg. Our theoretical light curves show excellent agreement with the optical y and infrared J, H, K 
light curves. The continuum UV 1455 A light curves observed with IUE are well reproduced. The turn-on and 
turn-off of supersoft X-ray observed with ROSAT are also explained simultaneously by our model. The WD 
mass is estimated, from the light curve fitting, to be M WD ps 1.15 M Q for V1500 Cyg, M WD rj 0.95 M Q for 
V1668 Cyg, and Mwd ~ 0.95 - 1 .05 M Q for V1974 Cyg, together with the appropriate chemical compositions 
of the ejecta. 

Subject headings: novae, cataclysmic variables — stars: individual (V1500 Cygni, V1668 Cygni, V1974 
Cygni) — white dwarfs — X-rays: stars 



1. INTRODUCTION 

It has been widely accepted that a classical nova is a result 
of thermonuclear runaway on a mas s-accreting w hite dwarf 
(WD) in a close binary system (e.g., Warnen fl995l for a re- 
view). Theoretical studies have elucidated full cycles of nova 
outbursts, i.e., from t he mass accretion stage to the end of 
a nova outburst (e.g.. iPrialnik & Kovetz I \l 9951) and devel- 
opme nts of light curves for various speed class es of novae 
(e.g.. iKato & Hachisu l[l994l iKato 111991 11997b . However, 
detailed studies of individual objects such as light curve anal- 
ysis have not been fully done yet except for the recurrent no- 
vae (e.g.. lHachisu & Kato"ll2000albl BOOlalbl l2003al l2006al 
lHachisu et al. 1120001 120031) . In order to further develop a 
quantitative study of nova light curves, a new approach is re- 
quired, from which we can derive the WD mass, ejecta mass, 
and duration of the nova outburst. 

On the other hand, some characteristic features and com- 
mon properties of nova light curves have been suggested, 
the reasons for which we do not yet fully understand. For 
example, infrared (IR) light curves of novae often show a 
decline law of F\ oc f~ Q , i.e., the flux at the wavelength 
A dec ays proportionally to a power of time from the out- 
burst. lEnnis et al. 1 (11977) suggested that near IR JHK fluxes 
of VI 500 Cyg (Nova Cygni 1975) are proportional to t~ 2 
in the early stages of nova explosions, while at later times 
they decay more steeply as f" 3 (see also. lKawara et al. "Il976t 
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Galla gher & Nev 1119761) . IWoodward et"al~l (119971) found that 
near-infrared light curves of VI 974 Cyg (Nova Cygni 1992) 
showed a p ower law of a ~ 1 . 5 in th e early 100 days. 

Recently, Hachisu & Kato (2005) explained these power- 
law declines as free-free emission from an optically thin 
plasma outside the photosphere. Here we further extend their 
approach and develop a light curve model that is widely appli- 
cable to optical and IR light curves of various types of novae. 
Section|2Jdescribes our basic idea and methods as well as our 
light curve models. We propose theoretical nova light curves 
having a universal decline law of Fx oc t~ a , regardless of the 
WD mass, chemical composition of the envelope, or observa- 
tional wavelength bands in optical and IR. Then we apply our 
method to three classical novae, i.e., V1500 Cyg in |5J V1668 
Cyg (Nova Cygni 1978) in H and VI 974 Cyg in §5\ Discus- 
sion follows in 5j6] and finally we summarize our results in 

2. MODELING OF NOVA LIGHT CURVES 

After the thermonuclear runaway sets in on a mass- 
accreting WD, its envelope expands greatly to R p h > 100 Rq 
and a large part of the envelope is ejected as a wind. Figure 
^illustrates a typical nova evolution, from the maximum ex- 
pansion of the photosphere to the end of hydrogen burning. 
At the very early expansion phase, the photosphere expands 
together with the ejecta. Then the photosphere lags behind 
the head of ejecta as its density decreases. Figure |2]depicts a 
schematic illustration of such a nova envelope, in which free- 
free emission from an optically thin plasma contributes to the 
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FIG. 1 . — Evolution of nova outbursts: (a) after a nova explosion sets in, 
the photosphere expands greatly up to > 100 Rq, and the companion star is 
engulfed deep inside the photosphere; (b) after the maximum expansion, the 
photospheric radius shrinks with time and free-free emission dominates the 
flux at relatively longer wavelengths; (c) a large part of the envelope matter is 
blown in the wind and the photosphere moves further inside; (d) the compan- 
ion eventually emerges from the WD photosphere and an accretion disk may 
appear or reestablished again; (e) the photosphere further shrinks to a size of 
< 0.1 Rq; (f) the optically thick wind stops; (g) hydrogen nuclear burning 
stops and the nova enters a cooling phase. Hard X-rays may originate from 
internal shocks between ejecta (or a bow shock between ejecta and the com- 
panion) from stage (a) to (/) as indicated by a dashed line. The ultraviolet 
(UV) flux dominates from stage (b) to (/). Then the supersoft X-ray flux 
replaces the UV flux from stage (f) to (g). 



optical and IR continuum fluxes. The photosphere eventually 
begins to shrink at or around the optical maximum. The nova 
envelope settles in a steady state after the optical maximum 
until the end of hydrogen shell burning (see Fig. 

2.1. Optically thick wind model 

The decay phase of novae can be w ell represe nted wit h a se- 
quence of steady-state solutions (e.g.. lKato & Hachisu 11994 
and references ther ein). Using the same met hod and numeri- 
cal techniques as in Kato & Hachisu ( 1994), we have calcu- 
lated theoretical models of nova outbursts. 
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FIG. 2. — A schematic configuration of our nova ejection model: A large 
part of the initial envelope mass is ejected by the winds, which are accelerated 
deep inside the photosphere. After the optical maximum, that is, after the 
maximum expansion of the photosphere, the photosphere begins to shrink 
whereas the ejecta are expanding. The optically thin layer emits free-free 
radiation at relatively longer wavelengths while blackbody radiation from the 
photosphere dominates at shorter wavelengths. 



We have solved a set of equations for the continuity, equa- 
tion of motion, radiative diffusion, and conservation of en- 
ergy, from the bottom of the hydrogen-rich envelope through 
the photosphere (see Fig. [2}, under the condition that the so- 
lution goes through a critical point of steady-state winds. The 
winds are accelerated deep inside the photosphere, so they are 
called "opticall y thick winds." We have used updated OPAL 
opacities (Iglesias & Rogers 1996). We simply assume that 
photons are emitted at the photosphere as a blackbody with 
a photospheric temperature of 7p n . We call this the "black- 
body light curve model" to distinguish these from the "free- 
free emission model," which will be introduced later. 

The wind mass l oss rate, M w j n H, is obtained as an e igenvalue 
of the equations llKato & Hachisu 1 119941 lHachisu & Kato I 
2001b) if the WD mass (MwdX envelope mass (AM env ), and 
chemical composition (X,F,XcNOj^Ne,Z) are given. Here X 
is the hydrogen content, Y is the helium content, Xqno is the 
abundance of carbon, nitrogen, and oxygen, X^ e is the neon 
content, and Z = 0.02 is the heavy element (heavier than he- 
lium) content, in which carbon, nitrogen, oxygen, and neon 
are also included with the solar composition ratios. The nu- 
clear burning rate M nuc , photospheric radius 7? pn , photospheric 
temperature r pn , and photospheric luminosity L pn as well as 
the photospheric wind velocity v p h are also calculated as a 
function of the WD mass (Mwd), chemical composition of the 
envelope (X, Y, Xqno, X$ e , Z), and envelope mass (AM env ), 
i.e., 

M wind = f x ( AM env ,X,Y, Xqno ,X Ne ,Z, Mwd), ( 1 ) 
M nuc = f 2 ( AM env ,X,Y, X CN0 ,Xx e ,Z, Mwd) , (2) 

^ph =/3(AM en v,X,y,XcNO,^ r Ne,Z,MwD), (3) 

r ph =/ 4 (AM env ,X,F,XcNO,^Ne,Z,M W D), (4) 

Lph = fs (AM env ,X,Y, XqnO , ^Ne , Z, M W d) , (5) 

Vph = / 6 (AM env ,X,F,X C NO,XNe,Z,M WD ). (6) 

The physical proper ties of wind solutions have been exten- 
sively discussed in Kato & Hachisu ( 1994) and some exam- 
ples of the wind solutions have been published in our previous 
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TABLE 1 

Chemical abundance of classical novae 



\JU Id* l 


H 


CNO 


Ne 


Na-Fe 


V382 Vel 1999 


0.47 


0.0018 


0.0099 


0.0069 


V382 Vel 1999 


0.66 


0.043 


0.027 


0.0030 


CPCru 1996 


0.47 


0.18 


0.047 


0.0026 


V723 Cas 1995 


0.52 


0.064 


0.052 


0.042 


V1425 Aql 1995 


0.51 


0.22 


0.0046 


0.0019 


V705 Cas 1993 #2 


0.57 


0.25 




0.0009 


V4169Sgr 1992 #2 


0.41 


0.033 






VI 974 Cyg 1992 


0.55 


0.12 


0.06 




V1974Cyg 1992 


0.19 


0.375 


0.11 


0.0051 


VI 974 Cyg 1992 


0.30 


0.14 


0.037 


0.075 


V351 Pup 1991 


0.37 


0.32 


0.11 




V838 Her 1991 


0.60 


0.028 


0.056 




NovaLMC 1990 #1 


0.18 


0.75 


0.026 


0.014 


V443 Set 1989 


0.49 


0.060 


0.00014 


0.0017 


V977 Sco 1989 


0.51 


0.072 


0.26 


0.0027 


V2214 Oph 1988 


0.34 


0.37 


0.017 


0.015 


QV Vul 1987 


0.68 


0.051 


0.00099 


0.00096 


V827 Her 1987 


0.36 


0.34 


0.00066 


0.0021 


V842 Cen 1986 


0.41 


0.36 


0.00090 


0.0038 


V842 Cen 1986 


0.58 


0.049 




0.0014 


QU Vul 1984 #2 


0.638 


0.034 


0.034 


0.005 


QU Vul 1984 #2 


0.36 


0.26 


0.18 


0.0014 


QU Vul 1984 #2 


0.33 


0.25 


0.086 


0.063 


QU Vul 1984 #2 


0.30 


0.06 


0.040 


0.0049 


PWVul 1984 #1 


0.62 


0.13 


0.001 


0.0027 


PWVul 1984 #1 


0.47 


0.30 


0.0040 


0.0048 


PWVul 1984 #1 


0.69 


0.066 


0.00066 




PWVul 1984 #1 


0.49 


0.28 


0.0019 




GQ Mus 1983 


0.37 


0.24 


0.0023 


0.0039 


GQ Mus 1983 


0.27 


0.40 


0.0034 


0.023 


GQ Mus 1983 


0.43 


0.19 






VI 370 Aql 1982 


0.044 


0.28 


0.56 


0.017 


V1370Aql 1982 


0.053 


0.23 


0.52 


0.11 


V693 CrA 1981 


0.40 


0.14 


0.23 




V693 CrA 1981 


0.16 


0.36 


0.26 


0.030 


V693 CrA 1981 


0.29 


0.25 


0.17 


0.016 


VI 668 Cyg 1978 


0.45 


0.33 






VI 668 Cyg 1978 


0.45 


0.32 


0.0068 




VI 500 Cyg 1975 


0.57 


0.149 


0.0099 




VI 500 Cyg 1975 


0.49 


0.275 


0.023 




HR Del 1967 


0.45 


0.074 


0.0030 




DQ Her 1935 


0.27 


0.57 






DQ Her 1935 


0.34 


0.56 






RR Pic 1925 


0.53 


0.032 


0.011 




TAur 1891 


0.47 


0.13 







reference 



Aueusto & Diaz (2003) 
Shore et al. TSoTB) 
l.vkcelal. (20031 
Iiuma (2006) 
LvkeetaT (2001) 
Aikhinmq ,-l <:nn*l i 
Scott etal. (W5) 
Vanlandinsham et al. (2005) 

Havward et^l^!996) 
Saizaretal (1996) 
Vanlandineham et al. 
Vanlandineham et al. , 
TndT^a^taT^^^ 
Andrea et al. 
Andrea et al. 
Andrea et al. 
Andrea et al. 
Andrea et al. 



(1997) 
(1999) 



' (19W) 
' (1994) 
S*J) 



(1994) 
(1994) 

i deFreitas Pacheco et al. ( 1989) 
Schwarz (2002) 

Andrea etaF (1994) 
Saizaretal" 199*2) 
Schwarz etalMTO7) 
Andrea etal. (1994) 
Saizaretal" 199*7) 
Andreae & Drechsel (19901 
Morisset & Pequienot (1996) 
Hassallet al. (1990) 
Andreae & Dr echsel (1990) 
Andrea etal. (1994) 
SniidersetaT (1987) 
Vanlandinsham et al. (1997) 
Andrea etal. (1994) 
Williams et al. (1985) 
Andrea et aTT 19§W 
Stickland etal. (1981) 
Lance etal. (1985^ 
Ferland & Shields (1978) 
Tvlenda (1978) 
Petitiean etaL (1990) 
Williams etaT (1978) 
Williams & Galla gto (1979) 
Gallagher et al. (1980a) 



TABLE 2 

Chemical composition of the present models 



novae case 


X 


^CNO 


^Ne 


Z a 


mixing b 


comments 


CO nova l d 


0.35 


0.50 


0.0 


0.02 


100% 


DQ Her 


CO nova 2 d 


0.35 


0.30 


0.0 


0.02 




GQ Mus 


CO nova 3 


0.45 


0.35 


0.0 


0.02 


55% 


VI 668 Cyg 


CO nova 4 


0.55 


0.20 


0.0 


0.02 


25% 


PWVul 


Ne nova l d 


0.35 


0.20 


0.10 


0.02 




V351 Pup, V1974 Cyg 


Ne nova 2 


0.55 


0.10 


0.03 


0.02 




V1500Cyg, V1974 Cyg 


Ne nova 3 


0.65 


0.03 


0.03 


0.02 




QU Vul 


Solar 


0.70 


0.0 


0.0 


0.02 


0% 





a carbon, nitrogen, oxygen, and neon are also included in Z = 0.02 with the same ratio as the 
solar abundance 

k ratio of mixing between the core material and the accreted matter with the solar abundances 
c chemical composition in the left columns is adopted for each nova listed below, although DQ 

Her, GQ Mus, PW Vul, V35 1 Pup, and QU Vul are discussed in separate papers 
^these three cases, CO nova 1, CO nova 2, and Ne nova 1, are hardly different from each other 

in their free-free light curves during the wind phase. Therefore, only the case of CO nova 2 is 

shown in this paper 
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papers fe.g.JHachisu & Katol2001albll2003bldl2004il2006a[ 
Hachi su et al. 1 1 996t 1 1 999alhl 1200(1 l2003t iKato ■ 1 9831 ll 997t 
1999). It should be noted that a large number of meshes, i.e., 
several thousand grids, are adopted when the photosphere ex- 
pands to R p h ~ 100 Rq. 

Using numeric tables of solutions Q-(|6|i with a linear in- 
terpolation between the adjacent envelope masses, we have 
calculated an evolutionary sequence by decreasing the enve- 
lope mass as follows: 

— AM env = Macc ~ Mwind ~ Mnuc , (7) 

at 

where M acc (< 10~ 9 Mq yr" 1 ) is the mass accretion rate onto 
the WD. We assume M ac? = in our calculation, because usu- 
ally M acc <C M w i n( j and M acc <C M nuc for classical novae. The 
envelope mass is decreased by the winds and nuclear burn- 
ing. A large amount of the envelope mass is lost mainly 
by the winds in the early phase of nova outbursts, where 

nuc- 

Optically thick winds stop after a large part of the enve- 
lope is blown in the wind (Fig. The envelope settles into 
a hydrostatic equilibrium where its mass is decreased by nu- 
clear burning, i.e., M w i n( j = in equation Q. Here we solve 
the equation of static balance instead of the equation of mo- 
tion. When the envelope mass decreases to below the min- 
imum mass for steady hydrogen-burning, hydrogen-burning 
begins to decay. The WD enters a cooling phase, in which the 
luminosity is supplied with heat flow from the ash of hydro- 
gen burning. 

In the optically thick wind model, a large part of the enve- 
lope is ejected continuously for a re latively long period (e.g., 
IKato & Hachisu 1119941 iKaroll 19971) . After the maximum ex- 
pansion, the photosphere shrinks gradually with the total lu- 
minosity (L p h) being almost constant. The photospheric tem- 
perature (r p h) increases with time because of L p h = 4ttR!L(tT' J,, 
The main emitting wavelength of radiation moves from opti- 
cal to UV. This causes the decrease in the optical luminosity 
and increase in the UV. Then the UV flux reaches a maxi- 
mum. Finally the supersoft X-ray flux increases after the UV 
flux decays. These timescales depend on the WD parameters 
such as the WD mass and chemical composition of the enve- 
lope (IKato II 19971) . Thus we can follow the developments of 
the optical, UV, and supersoft X-ray light curves by a single 
modeled sequence of the steady wind solutions. 

We have calculated various models by changing these pa- 
rameters. It should be noted here that the hydrogen content X 
and carbon, nitrogen, and oxygen content Xcno are important 
parameters because they are the main players in the CNO cy- 
cle but neon (X^ e ) is not because it is not involved in the CNO 
cycle. 

Table ^ summarizes the chemical composition of nova 
ejecta. Although each nova shows a different chemical com- 
position, we choose seven typical sets of the chemical com- 
positions as tabulated in Table[2] 

2.2. Duration of supersoft X-ray phase 

A luminous supersoft X-ray phase is expected only in a 
late stage of nova outbursts (see Fig. 0. Before the opti- 
cally thick winds stop, the photospheric temperature is not 
high enough to emit superso f t X-rays (T n h < 150,000 K; see, 
e.g.. IKato & Hachisu 1119941 lKato1fT997l) . This is because 
the winds are driven by the strong peak of OPAL opacity at 
T ~ 150,000 K. Moreover, a part of the supersoft X-ray flux 



TABLE 3 
Durations of the wind and 
hydrogen burning for co 

NOVAE 2 a 



WD mass 


^wind 


OH -burning 


(Mq) 


(days) 


(days) 


0.50 


2540 


7820 


0.55 


1960 


5760 


0.60 


1390 


4260 


0.65 


1040 


3220 


0.70 


858 


2560 


0.75 


598 


1740 


0.80 


496 


1370 


0.90 


319 


757 


1.00 


218 


459 


1 1 A 
1.1U 


1 J / 


lJU 


1.20 


95 


144 


1.30 


53 


66 


1.35 


29.3 


33.9 


1.37 


20.9 


24.1 


a chemical 


composition 


of the enve- 


lope is X = 


0.35, Xcno 


= 0.30, and 


Z = 0.02 








TABLE 4 




Durations of the wind and 


HYDROGEN BURNING FOR CO 




NOVAE 3 a 




WD mass 


^wind 


OH -burning 


(Mq) 


(days) 


(days) 


0.60 


1620 


5450 


0.70 


977 


3250 


0.80 


551 


1740 


0.90 


337 


941 


1.00 


234 


565 


1.10 


151 


303 


1.20 


100 


168 


1.30 


57 


77 


1.35 


33.2 


40.3 


1.37 


24.1 


28.0 


a chemical composition of the enve- 


lope is X = 


0.45, Zcno 


= 0.35, and 


Z = 0.02 








TABLE 5 




Durations of the wind and 


HYDROGEN BURNING FOR CO 




NOVAE 4 a 




WD mass 


^wind 


OH-burning 


(Mq) 


(days) 


(days) 


0.55 


3370 


11600 


0.60 


2392 


8528 


0.70 


1395 


5051 


0.80 


771 


2680 


0.90 


461 


1441 


1.00 


309 


851 


1.10 


196 


447 


1.20 


128 


241 


1.30 


70.0 


105 


1.35 


41.0 


52.7 


1.37 


30.0 


35.8 



a chemical composition of the enve- 
lope is X = 0.55, Xcno = 0.20, and 
Z = 0.02 
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may be self-absorbed by the wind itself. Here we roughly 
regard that a supersoft X-ray phase is detected only after the 
optically thick winds stop (X-ray turn-on). However, it should 
be noted that a recent X-ray observation of the 2006 out- 
burst of the recurrent nova RS Oph showed an earlier ap- 
pearance of the supersoft X-ray phase, the origin of which 
is not clear vet iHachisu & Katoll2006al iBode et al. 1120061 
Osbo rne et al. I2006albl) ~ 

When hydrogen shell-burning extinguishes, the supersoft 
X-ray flux drops sharply because the WD is quickly cooling 
down. This epoch corresponds to the turnoff of supersoft X- 
ray. In the present paper, we call stages (a) - (/) "the wind 
phase" because the optically thick winds blow during these 
stages, and stages (/) - (g) "the hydrogen-burning phase," be- 
cause the evolution is governed only by nuclear burning. The 
supersoft X-ray phase corresponds to the hydrogen-burning 
phase. 

We have calculated a total of 72 nova evolutionary se- 
quences. Tables |3]-|8] and Figure [3] show the duration of 
the wind phase and the epoch when hydrogen burning ends. 
The evolutionary speed of a nova depends on the WD mass 
and the chemical composition of the envelope. These two 
durations of the wind phase and of the hydrogen burning 
phase depend very weakly on the initial envelope mass, 
AM env ,o- Typical masses of AM env o are a few to several 
times 10~ 5 M Q and typical wind mass-loss rates are as large 
as M w i nc j ~ 10~ 4 - 10~ 3 Mm vr" 1 at a v ery early phase in 
our model (e.g., iKato & Hac hisu 1 U991 . The difference in 
AM env makes only a small difference in the durations, i.e., 
Af ~ AM enVi o/M w ind ~ 0.01-0.1 yr, at most. 

Supersoft X-ray phases were detected for several classical 
novae, whic h are a mani festation of hydrogen shell burning on 
a WD (e.g.. lOrioll2004t for recent summary). A full duration 
of the supersoft X-ray phase was obtained by Krautte ret al. I 
( 1996) for V1974 Cyg 1992, that is, about 250 days and 600 
days after the outburst for the rise (turn-on) and fall (turnoff) 
times, respectively. These two epochs are plotted in Fig- 
ure |5Jl (large open squares). Comparing our model with the 
observation, we are able to determine the WD mass to be 
~ 1.05-1.10 Mq, for the chemical composition of X = 0.55, 
X CN0 = 0.10, X Ne = 0.03, and Z = 0.02 (case Ne 2 in Table 
Et. This will be show n in detail in §0 For GQ Mus 1983, 
Shanl evet al. I i ll 9951) detected only the fall of supersoft X- 
ray phase nearly a decade after the outburst. We have roughly 
determined the WD mass of ~ 0.65 M Q for the chemical com- 
position of X = 0.35, Xcno = 0.30, and Z = 0.02 (case CO 2 
in Table |3, as shown in Figure [3^. Thus, the turn-on/turnoff 
time of supersoft X-ray is a good indicator of the WD mass. 

2.3. Nova light curves for free-free emission 

The blackbody light curves do not well reproduce the ob- 
served visual light curv es of novae as already discussed by 
lHachisu & Kato I (12005). Instead these authors have made 
light curves for free-free emission from optically thin ejecta 
outside the photosphere that can be reasonably fitted with the 
optical light curve of V1974 Cyg 1992. Thus, the above au- 
thors conclude that, except a very early phase of the outburst, 
optical fluxes of novae are dominated by free-free emission of 
the optically thin ejecta as illustrated in Figure[2] 

The flux of free-free emission (thermal bremsstrahlung) is 



jvdQdVdtdv - 



_ 16 /7T 

= T V6 



1/2 e 6 z 2 



1/2 



8 exp 



hv 

We 



TABLE 6 
Durations of the wind and 
hydrogen burning for neon 

NOVAE 2 a 



WD mass 


^wind 


^H-burning 


(Mq) 


(days) 


(days) 


0.55 


4010 


13600 


0.60 


2890 


10100 


0.70 


1680 


5970 


0.80 


914 


3150 


0.90 


533 


1690 


1.00 


356 


996 


1.05 


280 


720 


1.08 


248 


599 


1.10 


222 


510 


1.15 


182 


382 


1.20 


145 


280 


1.30 


80 


121 


1.35 


46.0 


60.5 


1.37 


32.5 


39.7 


a chemical composition of the enve- 


lope is X = 0.55, Xcno = 


0.10, X Nc = 


0.03, and Z = 


:0.02 






TABLE 7 




Durations of the wind and 


HYDROGEN BURNING FOR NEON 




NOVAE 3 a 




WD mass 


^wind 


^H-burning 


(Mq) 


(days) 


(days) 


0.55 


6180 


22800 


0.60 


4530 


16600 


0.70 


2550 


9780 


0.80 


1360 


5140 


0.90 


778 


2730 


1.00 


504 


1600 


1.10 


305 


817 


1.20 


198 


427 


1.30 


106 


177 


1.35 


59.9 


83.9 


1.37 


43.0 


55.3 



a chemical composition of the enve- 
lope is X = 0.65, Zcno = 0.03, X Ke = 
0.03, and Z = 0.02 



TABLE 8 

Durations of the wind and 
hydrogen burning for the 
solar abundance 



N e NidVldVdtdv, 
(8) 



WD mass 


^wind 


*H— burning 


(Mq) 


(days) 


(days) 


0.60 


6570 


25600 


0.65 


4810 


18700 


0.70 


3760 


14800 


0.75 


2540 


9940 


0.80 


1960 


7860 


0.90 


1100 


4070 


1.00 


712 


2350 


1.10 


419 


1190 


1.20 


263 


612 


1.30 


140 


247 


1.35 


78.1 


113 


1.37 


50.9 


72.9 
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(a)' SSS Phase' (X=0.35, CN0=0.30, Z=0.02) 



A H-buming ends 
O wind stops 




X— ray turn— off 
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white dwarf mass (M e ) 
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(b) SSS Phase (X=0.45, CNO=0.35, Z=0.02) 
A H-buming ends 
O wind stops 

X— ray turn— off 




X— ray turn-on 
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white dwarf mass (M ) 
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(c) SSS Phase (X=0.55, CNO=0.20, Z=0.02) 
A H-buming ends 
O wind stops 

X— ray turn— off 




X— ray turn— on 



0.6 0.8 1 1.2 

white dwarf mass (M e ) 



1.4 



(d) SSS Phase (X=0.55, CNO=0.10, Ne=0.03, Z=0.02) 
A H-buming ends 
O wind stops 

X— ray turn— off 
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white dwarf mass (M ) 
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(e) SSS Phase (X=0.65. CNO=0.03, Ne=0.03, Z=0.02]i 
A H-buming ends 
O wind stops 

X— ray turn— off 



X-ray turn— o 
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white dwarf mass (M e ) 
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(f) SSS Phase (X=0.70, Z=0.02) 



A H-buming ends 
O wind stops 
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—ray turn— off 
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white dwarf mass (M ) 



FIG. 3. — The turn-on (circles) and turnoff (triangles) times of supersoft X-ray are plotted against the WD mass. X-ray turn-on corresponds to the epoch 
when optically thick winds stop while X-ray turnoff corresponds to the epoch when hydrogen shell-burning ends. A total of 72 cases are plotted for six different 
chemical compositions: (a) case CO 2 in TableEJ a large square indicates the epoch of supersoft X-ray turnoff for GQ Mus 1983 i Shanlev et al. 1995); (b) case 
CO 3; (c) case CO 4; (d) case Ne 2, two large squares indicate the epochs of supersoft X-ray turn-on and turnoff for V1974 Cyg 1992 i &a utter^talTl9 96 ): (e) 
case Ne 3; (/) solar composition. 



where j v is the emissivity at the frequency v, Q, is the solid 
angle, V is the volume, t is the time, e is the electron charge, 
Z is the ion charge in units of e, c is the speed of light, m e 
the electron mass, k is the Boltzmann constant, T e the elec- 
tron temperature, g is the Gaunt factor, h the Planck constant, 
and and N , are the number densities of electrons and ions 
lAllenll98lL p.1031. 

The electron temperatures of nova ejecta were suggested to 
be around T e ~ 10 4 K and almo st constant during the nova out- 
bursts fe.g.. lEnnis et al. 119771 for V1500 Cyg). If we further 



assume that the ionization degree of ejecta is constant during 
the outburst, we have a flux of 

F\ oc / N e N t dV oc / ^^-?dr oc (9) 

for free-free emission of the optically thin ejecta during the 
optically thick wind phase, where Fx is the flux at the wave- 
length A, N e oc p w i n( j and Ni oc p w ; n d. Here, we assume v w i n( j = 
v p h and use the relation of continuity, p wind =M w ind/47rr 2 v w i nd , 
where p w j n d and v w ; n( i are the density and velocity of the wind, 
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(f) Solar abundance 
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FIG. 4. — Calculated light curves for free-free emission during the optically thick wind phase, (a) case CO 2 in Tabled (b) case CO 3, (c) case CO 4, (d) case 
Ne 2, (e) case Ne 3, (/) solar composition. Each panel shows light curves for different WD masses, i.e., Mwd = 1.3, 1.2, 1.1, 1.0, 0.9, 0.8, 0.7, and 0.6 Mq from 
left to right. The apparent magnitude constant is set to be = 2.5. 



respectively. We have obtained the absolute magnitude for 
free-free emission by 

(10) 

or the apparent magnitude by 

»A = "2.5 log ( f- d , V ( Vph , V 2 ( ^-\\c x , 
g VlO- 4 M Q yr-V Vl000kms-V \Rq J 

(11) 

where M\ (m\) is the absolute (apparent) magnitude for free- 



free emission at the wavelength A, and C\ (c\) is a constant. 
Subtracting equation dlOl from equation il It . we have 

c x = (m-M) x +Cx, (12) 

where (m-M)\ is the apparent distance modulus. In the 
present paper, we call c\ the apparent magnitude constant. 
We cannot uniquely specify the constant, C\ (or c\), because 
radiative transfer is not calculated outside the photosphere. 
Instead we choose the constant to fit the light curve. 

Figure |4] shows the free-free light curves calculated from 
equation with a fixed value of c\ = 2.5 for the six cases 
in Table |2 Here we plot the light curves only during the op- 
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FIG. 5. — The apparent magnitude constant, c\, in eq. 1111 is plotted against the apparent distance modulus, (m-M)x, for three novae, V1500 Cy g (O pen 
circle), V1668 Cyg (filled circ le), and V1974 Cyg (open square). Each panel shows (a) y (or V), (b) J, (c) H, and (d) K bands. Dash-dotted line: eq. 1121 for 
V1668 Cyg. Dashed line: eq. JlSj both for V1500 Cyg and V1974 Cyg. 



tically thick wind phase. The wind stops at the lower edge of 
each line. It is clear that the heavier the WD mass is the faster 
the development of a nova light curve is. Thus we conclude 
that the nova speed class is in principal closely related to the 
WD mass even when free-free emission dominates the nova 
optical fluxes. 

Using the individual fitting results described below in |3] 
13 and|5| we plot the apparent magnitude constant, c\, against 
the apparent distance modulus, (m—M)\, each for the y (or V), 
J, H, and K bands in Figure |5] If these three novae are just 
on a straight line with a slope of unity, we may conclude that 
the constant C\ for the absolute magnitude is universal among 
various novae. However, C\ for V1668 Cyg is about 1 .0- 1 .4 
mag dimmer than that for V1500 Cyg and V1974 Cyg. Each 
nova probably has a different C\, although we cannot draw 
any firm statements from only three novae. 



The distance modulus is calculated from 
(m—M)\ = 5 + 5logd+A\, 



(13) 



where d is the distance to the object. Absorption law for 
each band is given by A y = Ay = 3.1E(B -V) for y and V, 
Aj = O.SIE(B-V) for J, A H = 0.54E(B-V) for H, A K = 
0.35E(B-V) for A:, and/ U = OASE(B-V) for L band (e.g., 
IRieke & Lebofskvll985h . 

2.4. Template light curve of classical novae 

We have found that the free-free light curves calculated for 
various WD masses are quite homologous and overlap with 
each other as shown in Figure [6] Here we shift each light 
curve up and down and back and forth to overlap it with the 
1.0 Mq WD model, which is fixed in the figure. Table|5]lists 
the horizontal shifts in the m\ -logf plane that correspond to 
a time scaling factor of the light curve. 
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FIG. 6. — Same as those for Fig. |4] but each light curve is shifted vertically and horizontally to show a universal decline law. We can see a break on the 
free-free emission light curves near 100 days after outburst (200 days after outburst for the solar composition). We define this time as ?break- 



Figure Q shows that the 1 .0 M Q WD models with differ- 
ent chemical compositions are also homologous and can be 
overlapped with each other in the m\-\ogt plane except for 
the solar composition. Thus we may conclude that all these 
five light curves are homologous independently of the WD 
mass or the chemical composition. Only the solar composi- 
tion model shows a bit smaller decline in the m\ -logf plane. 
We call this property "a universal decline law of classical no- 
vae." Once a template light curve is given, we can specify a 
nova light curve using only one parameter and we choose the 
time at the "knee" as such a parameter. The template light 
curve has a prominent knee at logf w 2 (f ~ 100 days after 



outburst) in Figure 05, which is the only parameter that char- 
acterizes the timescale of nova light curves. This break time, 
fbreaki is listed in Table fTol for each set of the WD mass and 
chemical composition. 

It should be noted that this universal decline rate, a = 
1.7-1 .75 (here, F\ oc t~ a , t is the time in days after outburst), 
yields a simple relation between ?2 and t$, where ?2 is the time 
in days during which the star decays by 2 mag from the opti- 
cal maximum and is the time in days in 3 mag decay from 
the optical maximum. Since the flux obeys Fx oc t~ a before 
the break, we have 

m\ = -2.5 log Fa + const. = 2.5alogf + const. (14) 
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TABLE 9 

Time scaling factor of the nova light curves* 



WD mass 


CO 2 


CO 3 


CO 4 


Ne2 


Ne 3 


solar 


(Mq) 














0.60 


0.65 


0.70 


0.76 


0.79 


0.82 


0.87 


0.70 


0.47 


0.50 


0.57 


0.57 


0.60 


0.60 


0.80 


0.30 


0.32 


0.33 


0.33 


0.36 


0.34 


0.90 


0.12 


0.14 


0.13 


0.14 


0.14 


0.14 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


1.10 


-0.20 


-0.17 


-0.21 


-0.17 


-0.17 


-0.22 


1.20 


-0.40 


-0.35 


-0.39 


-0.34 


-0.34 


-0.41 


1.30 


-0.64 


-0.62 


-0.67 


-0.61 


-0.62 


-0.66 


a time scaling factor is in 


a logarith 


mic form of log £ 







TABLE 10 

Time at the break for the free-free emission light 

CURVES 



WD mass 


CO 2 


CO 3 


CO 4 


Ne2 


Ne3 


solar 


(Mq) 


(days) 


(days) 


(days) 


(days) 


(days) 


(days) 


0.60 


390 


436 


616 


812 


1120 


1510 


0.70 


257 


275 


398 


490 


676 


832 


0.80 


174 


182 


229 


282 


389 


457 


0.90 


114 


120 


145 


182 


234 


288 


1.00 


87 


87 


107 


132 


169 


209 


1.10 


55 


59 


66 


89 


115 


126 


1.20 


35 


39 


44 


60 


78 


81 


1.30 


20 


21 


23 


32 


41 


46 



The decay parameters of light curves, f 2 and f 3 , are related to 
the parameter a as 

1 = 2.5a [log(f 3 + Afo)-log(f 2 + Afo)] , (15) 

where Afo is the time in days from the outburst to the optical 
maximum. If a = 1.75, we have 



f 3 = 1.69*2+0.69 Af 



(16) 



Since the rise time Afo is usually short, from a few days (fast 
novae) to several days (moderately fast novae), the relation 
between f 2 and f 3 is very consistent with an (observational) 
empirical relation 

t 3 = (1. 68 ± 0.08) t 2 + (1.9 ±1.5) days, forf 3 < 80 days (17) 

or 

r 3 = (1.68±0.04)f 2 +(2.3±1.6)days, for f 3 > 80 days, (18) 
given by Cap accioli et al. 

2.5. UV 1455 A light curve 

Cassatella et al. I (|2()02) adopted two UV continuum bands 
to describe the UV fluxes of classical novae based on the IUE 
observation. One is the UV 1455 A band with a 20 A width 
(centered on 1455 A ) and the other is the UV 2885 A band 
with a 20 A width (centered on 2885 A ). Both of the bands are 
selected to avoid prominent emission or absorption lines in the 
UV spectra. In ourblackbody light curve model, the UV 1455 
A flux reaches its maximum at a photospheric temperature of 
~ 27,000 K while the UV 2885 A band reaches its maximum 
at a lower photospheric temperature of ~ 13,000 K. 

In this paper, we adopt only the UV 1455 A band be- 
cause our blackbody light curve model for the UV 1455 A 
band nicely follows the UV flux near the temperature of 
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FIG. 7. — (a) Theoretical free-free light curves for the 1.0 Mq WDs with 
six different chemical compositions, (b) The above six light curves are shifted 
vertically and horizontally to be overlapped with each other. Only the case of 
solar abundance (dashed tine labeled by "solar") is somewhat different from 
the other five cases. We can see a break on the free-free emission light curves 
near 100 days after outburst from a slope of ~ t 75 to ~ 5 . 



~ 27,000 K. This is confirmed partly by the theoretical nova 
spectra calculated by Hauschi ldt et al. I j!995l) . They showed 
non-LTE spectra ranging from optical to UV in their Figure 3. 
Their UV 1455 A band region is well fitted with a blackbody 
spectrum of T e ff = 25,000 K, around which our UV 1455 A 
flux reaches its maximum. This suggests that the blackbody 
light curve model for the 1455 A band is a good approxima- 
tion to the nova UV flux near its flux maximum. This is ob- 
servationally confirmed by direct fittings with the IUE obser- 
vation for VI 668 Cyg in 5g|and for VI 974 Cyg in |5] 

On the other hand, our blackbody light curve model for the 
UV 2885 A band is not a good approximation to the nova 
UV flux at/near its maximum partly because our blackbody 
2885 A flu x is about two times larg er than the UV flux cal- 
culated bv lHauschildt et al. I \ 19951) at T ef{ = 10,000 K, near 
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FIG. 8. — The ultraviolet (UV) 1455A and free-free emission light curves during the wind phase for six different chemical compositions: (a) case CO 2, (b) 
case CO 3, (c) case CO 4, id) case Ne 2, (e) case Ne 3, and (/) the solar composition. Eight light curves for different WD masses are calculated in each panel, 
i.e., Mwd = 1.3, 1.2, 1.1, 1.0, 0.9, 0.8, 0.7, and 0.6 Mq from left to right. The horizontal timescale is linear (not logarithmic) but stretched or squeezed to overlap 
each other. This time scaling factor is the same as those in Tablel9l The magnitudes of free-free light curves are arbitrary shifted to overlap with each other. The 
vertical flux scale for the UV 1455A is linear and normalized by the maximum flux. 



which the 2885 A flux reaches its maximum. To summarize, 
avoiding strong emission/absorption lines, we use the 1455 A 
band as a UV evolution of a classical nova, which is acciden- 
tally well followed by the blackbody light curve model near 
its flux maximum. 

Our light curves of the UV 1455 A band are plotted together 
with the free-free emission light curves in Figure|8]for various 
WD masses and chemical compositions. Note that the hori- 



zontal axis is not logarithmic but linear in this figure. These 
are stretched or squeezed by the same factor tabulated in Ta- 
ble|9] Light curves almost overlap with each other. Therefore, 
we may conclude that the UV 1455 A light curve is also spec- 
ified by only one parameter, such as a time scaling factor or 

^break ■ 

2.6. Relations among various nova timescales 
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FIG. 9. — Four typical nova timescales are plotted against the WD mass for six different chemical compositions of the nova envelope. Each panel shows (a) 
fH-buming, the time when hydrogen burning ends, (b) f w i nt j , the time when optically thick winds stop, (c) A?ijv.fwhm, the UV 1455A duration defined by the full 
width at the half maximum, and id) fb re ak> the time of break. These timescales depend not only on the WD mass but also on the chemical composition. 



Figure [9] shows the various timescales that characterize 
nova outbursts: (a) fH-buming, the epoch when hydrogen shell- 
burning ends; (b) f w i nc j, the epoch when optically thick winds 
stop; and (c) Afrjv.FWHM, the UV 1455 A duration, which 
is properly defined by the full width at the half maximum 
(FWHM) of the UV 1455 A light curve. In addition to these 
three timescales, it shows (d) the time of break, fbreak, where 
the free-free light curve has a knee in the m\-\ogt plane. It is 
clear that these timescales depend not only on the WD mass 
but also on the chemical composition. 

Figure ^3 shows the first three timescales, fH-buming, fwind, 
and Afuv.FWHM against the last one, fbreak- Now these three 
former timescales are monotonic functions of fbreak, regard- 
less of the WD mass or the chemical composition. Therefore, 
once we determine fb ie ak from observations, we can predict the 
duration of a luminous supersoft X-ray phase of a nova, i.e., 
its turn-on (f w ind) and turnoff (?H-buming) times. We have added 
the corresponding epochs for three individual novae, VI 500 
Cyg, V1668 Cyg, and V1974 Cyg. 

2.7. Transition in classical nova light curves 

Figure ^2 shows a schematic light curve of the free-free 
emission. The early slope of F\ oc f~ L75 changes at fb rea k to 
the late slope of Fx oc t~ . This break occurs a bit before the 
wind stops. This change of slope is caused mainly by a sharp 
decrease in the wind mass-loss rate, M w ; nc j, together with an 
increase in the wind velocity, v p h. 



After the optically thick wind stops, the free-free emission 
light curve changes its decline rate because no additional mass 
is supplied to the ejecta. In such a case, the free-free flux can 
be roughly estimated as 

f M 2 
F x oc / N e NidV oc p 2 V oc -^|V oc R~ 3 oc f~ 3 , (19) 

(e.g.. IWoodward et al. l[l9 97), where p is the density, M e j is 
the ejecta mass (M e j is constant in time after the wind stops), 
R is the radius of the ejecta (V oc R 3 ), and t is the time after 
the outburst. Therefore, in the later phase, the free-free flux 
decays as f~ 3 . In the actual case, the transition from f~ 35 to f~ 3 
may occur just before the wind stops. 

3. V1500 CYG (NOVA CYGNI 1975) 

We will examine now individual novae. The first example is 
an extremely fast nova V1500 Cyg. It has probably the fastest 
and largest eruption among novae. It rose to a maximum of 
my = 1.85 o n 1975 August 31 f rom a preoutburst brightness 
of m v > 21 jYoung et al. 11197(3) . A distance of 1 .2 ± 0.2 kpc 
tLance et aL _[lg88) and an interstellar extinction of E(B— V) = 
0.5+0.05 JFerland 119771) suggest a peak absolute luminosity 
of My = -10.0±0.3, which is about 4 mag brighter than the 
Eddington luminosity for a 1.0 M WD (see, e.g., eq. (4) of 
lHachisu & Kato II2004L for the Eddington V magnitude). An 
extensive summary of the observati onal results and modelings 
can be found in the review bv lFerland et all 11986). 
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FIG. 10. — Three typical nova timescales of 'H-buming, 'wind, and 
A'uv.FWHM are plotted against 'b,- C ak . These three timescales are monotonic 

functions of 'b re ak> so we can predict the epochs of 'H-burnine ™d 'wind if 'break 

is observationally determined. Estimated values of 'break are indicated by 
arrows for V1500 Cyg, V1668 Cyg, and V1974 Cyg. Large open squares 
denote the corresponding timescales determined from the observations. 




log[days after outburst] 

FIG. 1 1 . — A schematic nova light curve for our free-free emission model. 
When strong emission lines such as [O III] emission lines begin to contribute 
to the nova light curve, its light curve is gradually lift up (thin solid) from our 
template nova light curve (thick solid). This occurs when the nova enters a 
nebular phase. We recommend observations with narrower band filters such 
as the Stromgren y band to accurately follow the continuum flux even in the 
nebular phase. 



3.1. Light curve fitting in the decay phase 

iGallagher & Nevl (U976I) obtained the magnitudes of the 
three broad optical V, R, and / bands and the eight infrared 
1.2, 1.6,2.2,3.6, 4.8,8.5, 10.6, and 12.5 /im bands during the 
50 days following the discovery. They estimated the outburst 
day to be August 28.9 UT from the data of angular expansion 
of the photosphere. They conclude that the spectrum energy 
distribution is approximately that of a blackbody during the 
first 3 days while it is close to F v = constant after the fourth 
day. This F v = constant spectra resemble those usually as- 
cribed to the free-free emission. 

Based on the i nfrared photometry from 1 to 20 /zm, 
lEnnis et al. I dl977[> also presented s i milar observational re- 
sults like those of Gallagher & Nev (1976). The nova spec- 
trum changed from a blackbody to a bremsstrahlung emission 
at day ~ 4-5, that is, from that of a Rayleigh- Jeans (F„ oc v 2 ) 
to that of a thermal bremsstrahlung emission (F v ~ constant). 
Thus we regard that the nova enters a phase, in which free-free 
emission dominates, about 5 days after the outburst. 

They also obtained the onset of outburst on 
JD 2,442,653. 0±0. 5 from an ana lysis of the phot ospheric 
expansion similar to that by IGallagher & Nev I lll976l) . 
Therefore we define the outburst day of VI 500 Cyg as 
JD 2,442,653.0 (f = 0) in our model. 

Figur e IT2l shows the y ma gnitude light curve observed 
by |Lockwood & Millis I (fl976T) together with the three near- 
IR J (1.2 um), H (1.6 urn), and K (2 . 2 um ) magn itudes 
(Gallagher & NeJ] Il976t iKawara et al l Il976t iHnnis et.al.l 
119771) . We have measured the time of "knee" on these four 
light curves and obtained an average value of f bre ak ~ 70 days 
after the outburst. This value indicates a WD mass of 1 . 15 M Q 
among 1.1, 1.15, and 1.2 M Q WDs (see Table [Tol. Here 
we assume chemical composition of X = 0.55, Xcno = 0.10, 
= 0.03, and Z = 0.02 (c ase Ne 2 in Table 0, which is 
close to the estimate given bv lLance et al. I 1H988) as tabulated 
in TableG] 

Our best-fit model closely follows all the observations of 
optical y and infrared J, H, and K bands during the period 
from several to ~ 100 days after the outburst as seen in Fig- 
ure[21 We have obtained c y = 4.25, cj = 1.67, ch = 1.77, and 
ck = 1 .32 to fit with the V1500 Cyg light curves as already 
introduced in Figure [5] These well-fitted results confirm that 
free-free emission dominates the optical and infrared contin- 
uum at least from several days to ~ 100 days after the out- 
burst. 

Figu reH2ll also shows the light curve model by Enni s et al. I 
( 119771) . They assumed an expanding shell with an initial 
thickness of H and the velocity of v. Then the time depen- 
dence of the flux is approximately given by 



F v (xr 2 (l+2c s t/H) 



(20) 



In early stages of the expansion, the flux will be proportional 
to f" 2 , while it is more close to ?~ 3 in later times. They as- 
sumed a temperature of 10,000 K that gives a sound velocity 
of c s ~ 10 km s" 1 . The transition at ~ 60 days between the two 
slopes derives H ~ 10 13 cm, which indicates a shell-ejection 
duration of a half day after the explosion with an expansion 
velocity of v ~ 2000 km s -1 . The resultant light curve is not so 
closely following the observational data compared with that 
of our best-fit model. 

The wide-band V magnitude observations are plotted in 
Figure ^] In general, the wide-band V magnitude suffers 
from contamination of strong emission lines especially in the 
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FIG. 12. — Calculated free-free light curves are plotted for the three WD mass models: 1.1 Mq (dash-dotted lin e), 1.15 Mq (thick solid line), and 1.2 Mq 
(medium thickness solid line). Our free-free light curves nicely follow the observation after day ~ 5 as suggested by Ennis et al. 1 1 1977I) . Comparing /break with 
the observation, we choose the 1.15 Mq model among the 1.1, 1.15, and 1.2 Mq models (see Table llOl . We also added light curves of F\ oc t~ 3 after the optically 
thick wind stops. Here we assume a chemical com position of X = 0.55, Xqno = 0.10, Xn c = 0.03, and Z = 0.02 (case Ne 2 in Table |2*1. Each pane l shows (a 
magnitude (open circle: Lockwood & Millis 1976). (b) J, (c) H, and (d) K magnitudes. The data of near infrared J, H, and K bands are taken from E nnis et 
1 1977, open circles), Kawara et al. 1 1976| open diamonds), andfGallaeher & Nev 1 1976, filled circles). Thin solid line in (d): a model light curve of eq 
proposed b v IEnnis eta!TTl^7^i ~'^ 



nebular phase. Comparing with Figure fT2h. we see that the 
V magnitude decays a little bit more slowly than the y mag- 
nitude does during 10 < t < 100 days. This is due to the ap- 
pearance of strong [O III] 4958.9 and 5006.9 A emission lines 
located just at the shorter edge of the V bandpass. We see 
some systematic differences ev en in the three different V mag - 
nitudes among o pen squares t Arkhipova & Zaitseva 1976), 
open diamonds (Pfau 1976). and open circles dTemnesti 
[T979T) . This comes from the difference of the wavelength sen- 
sitivity among these three V systems near the shorter edge 
of the V bandpass. Although the visual magnitudes are the 
richest in observational points as shown in Figure[0] the dif- 
ference from the y band (even from the V band) becomes 
significant in the later nebular phase. Therefore we use the 



medium-band y magnitude for the optical light curve fitting. 
We must be careful that the V magnitude light curves do not 
precisely follow the continuum flux but are highly contami- 
nated by strong emission lines especially in later phases. 

3.2. Distance to V1500 Cyg 

The di stance to VI 500 Cyg has been discussed by many 
authors. Young et al. (1976) estimate d the distance to be 
1.4 ± 0.1 kpc for E(B-V) = 0.45 jTomkin et al. II 19761) from 
the reddening-distance law toward the nova, which is shown 
in Figure^] The data align almost in a straight line. Another 
estimate comes from an empirical relation between the max- 
imum absolute magnitude and the rate of decline (MMRD). 
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FIG. 13. — Same as Fig. 1121 but for the visual and V observations. 
Our model light curves are plotted for the three WD masses of 1.15 (thick 
solid), 1.1 (dash-dotted), and 1.05 Mq (thin solid). The visual observa- 
tional data are taken from AAVSO (s mall dots) an d the V data are taken 
from TempestD (1523, open circles), Pfajy 1197ft! open diamonds), and 
Arkhinova & Zaits eva|<19761 open squares). Three epochs are indicated by 
arrows for our best-fit 1.15 Mq model. 
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FIG. 14. — The distance-reddening law in the directi on of VI 50 
Cyg, each star (open circles) of which is taken from lYoung et al. 1 119761) . 
Solid line: maximum magnitude versus rate of decline relation (labeled 
"MMRD"). Dotted line with r everse carets: upper limit to (m—M)v derived 
bv lAndo & Yamashita I H976T) . 



Using Schmidt-Kaler's iSchmidtll957l) empirical relation 

My )max =-1 1.5+ 2.5 log* 3 , (21) 

we obtain M v _ mm » —10 .11 together with t-} = 3.6 days (e.g. 
iDuerbeck & Wolf 119771) . Combining the maximum apparent 
magnitude of my. max = 1.85 and the absolute maximum mag- 
nitude of Afymax = —10.11, we obtain a simple relation of 

(m-M) v = 5 + 5\ogd+ 3 AE(B-V)= 11.96. (22) 

This distance-reddening relation is plotted in Figure^] which 
is labeled "MMRD." 

A firm upper limit to the appa rent distance modulus was 
obtained to be (m-M) v < 12.5 bv lAndo & Yamashita1(ll976l) 
from the Galactic rotational velocities of interstellar H and K 
absorption lines. This upper limit is also plotted in FigurefTH 
which is labeled "AY." 

The nebular expansion parallax method is a different way 
to estimate the distance. Becker & Duerbeck ( 1980) first im- 
aged an expanding nebular (0"25 yr" 1 ) of V1500 Cyg and 



estimated the distance to be 1350 pc together wit h an expan- 
sion v elocity of v exp = 1600 km s" 1 . However, Wade et al~l 
J1991I) resolved an expanding nebular and obtained a much 
lower expansion rate of 0."16 yr" 1 . Using a much sma ller ex- 
pansio n velocity of v exp = 1180 km s" 1 observed by ICohen I 
( 1985) , they estima ted the distance to be 1.56 kpc. Finally, 
Slavi n et al. 1 JT995T) obtained a more expanding image of the 
nebular (0"16 yr" 1 ) and determined the distance to be 1550 
pc assuming v exp = 1 180 km s" 1 . 

Therefore we may conclude that the distance is d s=s 1 .5 kpc 
assuming a color excess of E(B-V) as 0.45 as shown in Fig- 
ure El Using this distance of d as 1.5 kpc and color excess 
of E(B-V) as 0.45, we have estimated the apparent distance 
moduli of {m-M)\ for the y, 7, H, and K bands, as already 
introduced in Figure|5] 

3.3. The white dwarf mass 

The WD mass of 1 . 15 Mq for our b est-fit model is consis- 
tent w ith Mwd > 0.9 Mq obtained bv lHorne and Schneider I 
( 1989), based on a spectral line analysis. Such a massive WD 
is consider ed to be an o xygen-neon-magnesium (ONeMg) 
WD. Umeda et al. ( 1999) obtained an upper limit for the 
mass of a CO WD in a binary, Mco < 1.07 Mq, when it is 
just born. The WD mass decreases after every nova outburst 
because a part of the WD core is dredge d up by convection 
and b lown off during the outburst (e.g., Prialnik & Kovetzl 
119951) . Therefore, the 1.15 M WD is probably not a CO WD 
but an ONeMg WD. This conclusion is also consistent with 
the chemical composition of the ejecta as tabu lated in Tabled 
JFerland & Shields II 97ft I71n7e et al. 11 9881). 

3.4. Emergence of the secondary component 

The mass of the donor star (the se condary component) can 
be estimated from the orbital period. Semen iuk et al. I i 19951) 
determined an orbital period of P 01D = 0.1396 days (3.351 hr) 
from the orbital modulations with an amplitude of ~ 1 mag. 
Using Warner's (1995) empirical formula 
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we have M2 = 0.29 Mq. Then the separation is a = 1.28 Rq, 
the effective radius of the Roche lobe for the primary compo- 
nent (WD) is Rl = 0.64 Rq, and the effective radius of the 
secondary is R2 = 0.34 Rq. In our model, the companion 
emerges from the WD envelope when the photospheric radius 
of the WD shrinks to R p h ^1.0 Rq (see Fig. 0. This epoch 
is about day 50 in our best-fit model of Mwd = 1.15 Mq, but 
this time there is no transition of any kind in the optical light 
curve. 

4. V1668 CYG (NOVA CYGNI 1978) 

V1668 C yg wa s discovered on September 10.24 UT, 1978 
(M orrison 1978), 2 days before its optical maximum of 
mv,mnx = 6.04. Since no estimate of the outburst day has ever 
been given, we assume it to be JD 2,443,759.0. The out- 
burst of VI 668 Cyg was well observed with IUE in the ul- 
traviolet (UV) during a full period of the UV outburst (e.g., 
ICassatella et al. 1119791 IStickland et al II1981I). In a pioneer- 
ing work of a nova light-curve model. lKato I dl994l) presented 
light curve fitting of VI 668 Cyg based on the blackbody 
model. However, it is too simplified for classical novae. Here 
we present new light curves for the optical and IR observa- 
tions based on the free-free emission model. 
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FIG. 15.— UV 1455 A light curves for the best-fit 0.95 Mq WD model 
(thick solid line) together with the 1 .0 Mq (thin solid line labeled by 1 .0 Mq ) 
and 0.9 Mq (thin solid line labeled by 0.9 Mq) WD models. Here we as- 
sume a chemical composition of X = 0.45, Xqno = 0.35, and Z = 0.02 (case 
CO 3 in TableEJ. The UV 1455 A data (large open circles) are taken from 
ICassatella et al. lE00 2V A large part of the UV flux is absorbed by dust after 
the formation of an optically thin dust shell. 



4.1. UV 1455 A fluxes 

Figure [151 shows li ght curve fitting of ou r UV 1455 A flux 
with the IUE data (Cassatella et al. 1120021) . Here we adopt 
a chemical composition of X = 0.45, Xcno = 0-35, and Z = 
0.02 (case CO 3 in T abled after Stickland et al. I (119811) and 
lAndrea et al. i d! 994). Regarding the steep fall of the UV flux 
on JD 2,443,815 as caused by a dust absorption, we choose 
the best-fit model of M W d = 0.95 Mq among 0.9, 0.95, and 
1.0 Mq WDs. We have to simultaneously fit our model light 
curves with both the optical and UV 1455 A observations as 
shown in Figures [161 

The chemical composition of the ejecta suggests a CO WD. 
Our estimated WD mass ofMwD ~ 0.95 Mq is consistent with 
the upper mass li mit for a CO W D at its birth in a binary, i.e., 
Mco < 1 -07 M Q dUmedaetal. 119991) . 

4.2. Distance to V1668 Cyg 

The distance-redde ning law in the di rection of VI 668 Cyg 
was obtained by Slovak & Vogt ( 1979), although the number 
of stars is small and the data are scattered as shown in Figure 
IT71 They also obtained a reddening of E(B-V) = 0.38 from 
the interstellar f eature of K I (7699 A) and then a distance 
of d = 3.3 kpc. Duerbeck et al. ( 1980) criticized Slovak & 
Vogt's work and proposed a distance of d = 2.3 kpc from their 
newly obtained distance-reddening law and E(B-V) = 0.35, 
based on the same stars depicted in Figure [n] Assum- 
ing that the opt i cal maximum is the Eddington luminosity, 
IStickland et al. I ( 119811) estimated the distance to be d = 2.2 
kpc together with their E(B- V) = 0.40 from the 2200 A fea- 
ture. It is clear that these distance-reddening laws rely only on 
the three stars beyond 1 kpc in Figure ^1 Since no accurate 
law is d rawn from these little data, w e cannot judge the results 
of both iSlovak & VogTl ff979l) and iDiierheck et all ff980l) . 
As for the result of Stickland et al., we have no evidence that 
the maximum luminosity of V1668 Cyg is just the Eddington 
limit. 

The distance to the nova can be also estimated from 
the absolute magnitude at the optical maximum v ersus rate 
of decline (MMRD) relation. iKIare et al. I (119801) obtained 
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FIG. 16. — Free-free emission and UV 1455 A light curves for the best- 
fit 0.95 Mq WD model (thick solid line) together with the 1.0 Mq (medium 
thickness solid line) and 0.9 Mq (dash-dotted line) WD models. A straight 
line of Fx °c ' 3 is also added after the optically thick winds stop. We assume 
a chemical composition of X = 0.45, Xcno = 0.35, and Z = 0.02 (case CO 3 
in Table El. Visual observation (small dots) is taken from AAVSO. The UV 
1455 A data (large open circles) are taken from Cassatella et al. 1 200J). Two 
epochs are indicated by an arrow: one is the formation epoch of an optically 
thin dust shell and the other is the epoch when the optically thick winds stop. 
The v magnitude (o ven squares: [G allagher et al. I19 80bl) and V magnitude 
(small open circles: Mallama & Skillman 1979) observations are added. It 
is clear that the y magnitude departs from the visual magnitude after day 
~50. 
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FIG. 17. — The distance-reddening law in the direction of VI 66 8 Cyg, 
each star (open circles) of which is taken from Slovak & Voat 1 1979). Solid 
lines: the distance-reddening relation calculated from the UV 1455A flux 
fitting (labeled "UV 1455 A") and the maximum magnitude versus rate of 
decline relation (labeled "MMRD"). These three trends/lines merge into one 
at the point of E(B -V)m 0.4 and d fa 3.5 kpc. 



^v, max = -8 .0 ± 0.2 from Schmidt-Kale r's (iSchmidt 119571) re- 
lation (1211 together with f 3 = 24.3 days ( iMallama & Skillman 
1979). This gives a distance of d = 3.7 kpc together with 
mv,nm. = 6.04 and A v = 3AE(B-V) = 1.24. Using these val- 
ues, we get a distance-reddening relation for the MMRD rela- 
tion 

-8.0 = 6.04-5-51ogd-3.LE'(fi-V), (24) 

which is plotted in Figure[l7](labeled by "MMRD"). 

We also add another distance-reddening relation calculated 
from our UV 1455A flux fitting (labeled by "UV 1455 A "), 
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FIG. 18. — Free-free light curves for the best-fit 0.95 Mq WD model (thick solid line and a lift-up thin solid line) together with the 0.90 Mq (dash-dotted) and 
1.0 Mq (medium thickness solid) WD models. Each panel shows (a) y magnitude iGallaeher et al. 1980b), (b) V magnitude (open circles), data of which are 
taken from Mallama & Skillmanl ll979D . together with y magnitude (open square), (c) J, (d) H, (e) K (2.3/im), and (/) L (3.6/^m) magnitudes. The near infrared 
data of /, H, K, and L (open circles) are taken from Gehrz et al. 1 1980). We also added free-free light curves that follow the brightness after the dust formation 
for the three infrared bands in (d)-(f). Two epochs are indicated: one is the epoch when the optically thin dust shell formed and the other is the epoch when the 
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i.e., 

2.51og(5.8x 10 _I3 ) = 2.51og(1.6x 10- 12 )-51og(— ^— )-8.3£(B 

10 kpc 

(25) 

where F\ = 5.8 x 10~ 13 ergs cm" 2 s" 1 A" 1 is the calculated 
flux for the 1455A band at the distance of 10 kpc. These three 
trends/lines cross each other at the point of E(B-V) f» 0.4 and 
d 3.6 kpc. This reddening is very consistent with the red- 
dening of E(B-V) = 0.4 obtained bv lStickland et al. I (119811) 
from the 2200 A feature. H ere we assume A\ = &.3E(B-V) 
at A = 1455A JSeaton 119791) . 

In this paper, we adopt a distance of d = 3.6 kpc and a 
reddening of E(B-V) = 0.40. Then the optical maximum 
exceeds the Eddington luminosity by ~ 2 mag because of 
M v ,Edd = -5 .85 for the 0.95 M WD. 

4.3. Optical and infrared magnitudes 

Figure [T8l shows the y magnitude (Gallagher et al. 1ll980hl) 
as we ll as the J, H, K (2.3/i), andL (3.6/i) bands ( Gehrz ^Takl 
1980). In Figure 118b. w e also plot the V magnitude 
(Mall ama & Skillman 1119791) . The y magnitude decays faster 
than the V magnitude, because it does not include the [O III] 
emission lines that contribute to the V magnitude. The V 
magnitude begins to dev iate from the y magnitude after day 
~ 50. iKlare et al71 ([1980) found that the nebular phase started 
at least 53 days after the optical maximum. Our best-fit model 
of Mwd = 0.95 M Q nicely follows the y-band data. 

The four IR band magnitudes of J, H, K (2.3 fim), and L 
(3.6 /j,m) also nicely follow our free-free light curves until 
the formation of an optically thin dust shell. Dust formation 
started from day ~ 30 and each IR flux reached its maximum 
on day ~ 60. We cannot trace the formation of an optically 
thin dust shell on the y magnitude light curve. 

We obtain the fitting constants in equations il It and (I12> . 
i.e., c y = c v = 6.85, cj = 4.78, c H = 4.50, c K = 4.55, c L = 3.63 
as already shown in Figure [5] Each distance modulus of (m— 
M)\ is calculated using a distance of d = 3.6 kpc and color 
excess of E(B-V) = 0.40. 

4.4. Emergence of the secondary component 

iKaluznvl d!990l> reported that VI 668 Cyg is an eclipsing 
binary with an orbital period of 0.1384 days (3.32 hr). We 
estimate M2 = 0.29 M Q from equation i23\ . Then the sepa- 
ration is a = 1.2 Rq, the effective radius of the Roche lobe 
for the primary component (WD) is R* = 0.59 Rq, and the 
effective radius of the secondary is = 0-34 Rq. When the 
photospheric radius of the WD shrinks to /? p h ~ 1.0 Rq, the 
companion emerges from the WD envelope (see Fig. [Q. This 
epoch is ~ 100 days after the outburst in our best-fit model of 
Mwd = 0.95 Mq. There are no transition at this epoch in the 
optical and IR light curves. 

5. V1974 CYG (NOVA CYGNI 1992) 

V1974 Cyg was disc overed at my ~ 6.8 on February 19.07 
UT (JD 2.448.671 .57: ICollins lfT992h on the way to its op- 
tical maximum of my )max « 4.2 around February 22 (JD 
2,448,674.5). Since the outburst day is not accurately esti- 
mated, we assume it to be JD 2,448,670.0 in the present paper. 
This was the first nova ever observed with all the wavelengths 
from gamma-ray to radio, especially well observed with the 
X-ray satellite ROSAT and the UV satellite WE. ROSAT first 
detect ed the on and off of supersoft X-ray from the classical 
novae ( iKrautter et al. ll996tBaiman et al. I199H) . 
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FIG. 19. — Calculated free-free emission, UV 1455 A, and supersoft X-ray 
light curves are plotted for the best-fit 1 .05 Mq (thick solid line) WD model 
together with the 1 .0 Mq (medium thickness solid line) and 1 . 1 Mq (dash- 
dotted line) WD models. Fittings of the supersoft X-ray flux (or count rate) 
and of the UV 1455 A flux are shown in Figs. 1201 and |2*T1 separately. We 
assume a chemical composition of X = 0.55, Xcno = 0.10, Xn c = 0.03, and 
Z = 0.02 (case Ne 2 in Tablel2l. The visual observations (small dots) are taken 
from A AVSO. The V magnitudes (small open circle) are from Chochol et al. 
119931) . The supersoft X-ray data (open squares) are from IKrautter et aLI 
1 1996). The UV 1455 A data (large open circles) are from Cassatella et al. 
(2002). Two epochs, which are observationally suggested, are indicated by 
large upward arrows: one is the epoch when the optically thick winds stop 
and the other is the epoch when the hydrogen shell-burning ends. 



lHachisu & Kato I (|2005) presented model light curves for 
VI 974 Cyg in the visual, UV, and X-ray bands and ob- 
tained the best-fit parameters of M WD = 1.05 Mq, X = 0.46, 
and Xc no = 015 (for g i ven X N b = 0.05). Based on this 
model, Kato & Hachisu ( 2005) presented an optical light 
curve model for V1974 Cyg in the very early phase (i.e., the 
super- Eddington phase). Here, we present a light curve analy- 
sis for V1974 Cyg again, adding IR data to fitting. We assume 
a different chemical compositi on of X = 0.55, Xcno = 0.10 , 
X Ne = 0.03, and Z = 0.02 after IVanlandingham et all ( 120051) 
to examine how the estimated WD mass depends on the given 
chemical composition. This composition (case Ne 2 in Ta- 
ble|2jl is slightly different from that of Vanlandingham et al. 
but hardly affects the behavior of light curves because neon is 
not included in the CNO cycle and the difference in the CNO 
abundance is so small. 

5.1. Fitting with multiwavelength light curves 

We start our model fitting at shorter wavelengths. Figures 
U9H21l show the supersoft X-ray fluxes observed with ROSAT 
(Krautter et al. 1996) an d the UV 1455 A continuum fluxes 
observed with WE Jcassatella et al. 120021) . 

For a fixed abundance of the nova envelope, the only free 
parameter is the WD mass. We have calculated supersoft X- 
ray light curves for a wavelength window of 0.1 -2.4 keV 
with the blackbody model. The best-fit one is obtained for 
the WD mass of 1.05 M Q among 1.0, 1.05, and 1.1 M Q as 
shown in Figures ^] and |20| The supersoft X-ray emerged 
on day ~ 260 after the outburst and remained almost constant 
during the next ~ 300 days, and then decayed rapidly on day 
~ 600. The 1.05 Mq WD model shows a bit longer duration 
of its supersoft X-ray phase compared with the observation 
(see also Table|6|i. 

Our calculated X-ray fluxes in Figure|20|show that the more 
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FIG. 20. — Calculated X-ray fluxes (0.1 -2.4 keV) are plotted against 
time for various WD masses together with the ROSAT observation count 
rates (open and filled circles: taken from IKrautter et alntl996j). Open cir- 
cles: dominated by soft X-rays. Small filled circles: dominated by hard X- 
rays. Open squares and small up ward arrows: corrected X-ray fluxes and 
lower limits i Balman et al. 1998). The epoch of the optical maximum corre- 
sponds to JD 2,448,673.67, which is 2.67 days after the outburst. Thin solid 
lines: 1.0, 1.08 and 1.1 Mq WDs with the envelope composition of X = 0.55, 
Xcno = 0.10, X Ne = 0.03, and Z = 0.02 (case Ne 2 in Table^. Thick solid 
line: the best-fit model of 1.05 Mq WD. Same two epochs as in Fig. [19] arc 
indicated by arrows. 
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FIG. 21 . — UV 1455 A light curves for the best-fit 1.05 Mq (thick solid 
line) WD model together with the 1.0, 1.08 and 1.1 Mq (thin solid lines) 
WD models. We assume a chemical composition of X = 0.55, Xcno = 0.10, 
X Nc = 0.03, and Z = 0.02 (c ase Ne 2 in Table B. The UV 1455 A data (large 
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massive the WD, the shorter the duration of the supersoft X- 
ray phase. This is because a stronger gravity in a more mas- 
sive WD results in a smaller ignition mass. As a result, hy- 
drogen is exhausted in a shorter period (see, e.g.. iKato 119971 
for X-ray turnoff time). Therefore we have calculated X-ray 
and UV light curves for a 1 .08 M© WD, which is between the 
1.05 and 1.1 Mq WDs. Fitting becomes much better in the 
supersoft X-ray but worse in the UV as easily seen from Fig- 
ure |^ Therefore, we regard the 1.05 M Q WD model as the 
best-fit one at least for the assumed chemical composition. 

Comparin g our best-fit UV 1455 A model with the observa- 
tion (|£assatella et al. 1 120021) we obtain a distance of 1.7 kpc. 
Here w e again adopt the absorption law given by Seaton 
d!979h . A\ = S.3E(B-V) = 2.65, together with an extinction 



of E(B-V) = 0.32 estimated bv lChochol etaTI (fl997h . 

5.2. Free-free light curves for optical and infrared 

We then examine the V, J, H, and K light curves. Figures 
1X91 and I22h show the three different V magnitude observa- 
tions. Two of them cl e arly d epart from each other after day 
~ 80. ICho chol et al. (1993) explained that this systematic 
difference in V magnitude comes from the [O III] emission 
lines locating at the shorter edge of the V bandpass. A differ- 
ent instrumental system with a slightly different V bandpass 
allows variation of the V-magnitude. At the end of 1992 April 
(day ~ 75) strong [O III] 4958.9 and 5006.9 A emission lines 
appeared and the V brightness slightly rose, creating a small 
bump in the light curve. The discrepancy is more prominent 
in the visual light curve among observers in Figure Tl9\(small 
dot: taken from AAVSO). We find no y magnitude observa- 
tions of this object. Therefore, we here adopt the faintest one 
among various V magnitude data for the light curve fitting. 
Figure \19\ shows that our free-free light curve follows the V 
magnitudes until day ~ 80 and it begins to deviate after that. 

We have also fitted the J (1.25/im), H (1.6^m), and K 
(2.3/im) magnitudes in Figure |22] using the same models 
as in Figure 119 1 The 7, H, and K data are taken from 
I Wood ward et al. I l[l997), where they summarized their IR ob- 
servations and concluded that the V, J, H, and K light curves 
all showed an abrupt transition from a Fx oc f" 1 5 slope to a 
Fx oc f~ 3 slope at day ~ 170. Our free-free light curve well re- 
produces these observations until day ^80-100 and deviates 
after that. 

We have obtained the fitting constants in equation (^J, 
which are cy = 3.73, cj = 1.68, c# = 1.75, and ck = 1.50, as 
already shown in Figure[5] 

5.3. Distance to V1974 Cyg 

IChochol et al. I i ll 9971) extensively discussed the distance to 
VI 974 Cyg mainly based on the maximum magnitude versus 
rate of decline (MMRD) relations and concluded that the most 
p robable value is 1 8 kpc. 

lHachisu & Kato I J2005I) derived a distance of d = 1.7 kpc 
from a direct fit with the model UV 1455 A light curve. 
For a distance of 1 .7 - 1 .8 kpc, the peak lu minos ity is super- 
Eddingtonby ~ 1.8 mag. Kato & Hachisu (2005) proposed a 
mechanism of the super-Eddington luminosity and calculated 
a super-Eddington light curve for V1974 Cyg. From the fit 
with the UV 1455 A model they obtain a distance of 1.8 kpc 
and a 1 .7 mag super-Eddington lumi nosity. Their distance is a 
bit larger than the previous result by Hachis u"& Kato I J2005I) . 
because the model UV luminosity increases due to the super- 
Eddington effect. 

In this paper, we adopt a distance of d = 1.8 kpc, a color 
excess of E(B-V) = 0.32, and the same absorption laws as in 
V1500 Cyg. 

5.4. Dependence on the chemical composition 

Fi gures 1231 and l24l show light curve fittings for another set 
of the chemical composition (case Ne 1 in Ta ble 13, w hich is 
close to the composition obtained by H avward et al. I (119961) 
as shown in Table [T] We obtain the best-fit 0.95 M Q WD 
model among 1.0, 0.95, and 0.9 M WDs. In this 0.95 M Q 
WD model, f w ; n d = 269 days and fH-buming = 594 days. The 
direct fit with our UV 1455 A model indicates a distance of 
1 .6 kpc. 



20 



4) 
TJ 

'c 
a> 
o 

£ 
I 

> 



5 - 



10 - 



15 



-i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

(a) V1974 Cyg (Nova Cyg 1992) V mag 



Hachisu & Kato 




CNO=0.10, Ne=0.03 
. X=0.55, Z=0.02 

=(1.1). 1.05, (1.0) M G 
i i I i i_ 



wind stops 
(day ~280) 

I 




-« 5 I" 
TJ 

■ 'c 

E 
l 

-P10 h 



1 2 
log [days after outburst] 



15 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

(b) V1974 Cyg (Nova Cyg 1992) J mag 



wind stops 
(day ~280) 




CNO=0.10, Ne=0.03 
X=0.55, Z=0.02 



=(1.1), 1.05, (1.0) M 



1 2 
log[days after outburst] 



TJ 
3 

°c 
a> 
o 
E 
I 

x 



5 - 



10 - 



15 



-i 1 1 r 



— | 1 1 1 1 1 1 1 1 1 1 1 r 

(c) V1974 Cyg (Nova Cyg 1992) H mag 



wind stops 
(day ~280) 




CNO=0.10, Ne=0.03 
X=0.55, Z=0.02 
M WD =(1.1), 1.05, (1.0) M G 



- n, 5 " 



-10 - 



1 2 
log[days after outburst] 



15 



1 1 1 1 1 1 1 1 1 1 1 1 r 

(d) V1974 Cyg (Nova Cyg 1992) K mag 




o wind stops 
(day ~280) 

I 



CNO=0.10, Ne=0.03 
X=0.55, Z=0.02 
M WD =(1.1), 1.05, (1.0) M Q 



1 2 
log[days after outburst] 



FIG. 22. — Calculated free-free emission light curves for the best-fit 1.05 Mq (thick solid line) WD model together with the 1.1 Mq (dash-dotted) and 1.0 Mg 
(medium thickness solid) WD models. We assume a chemical composition of X = 0.55, Xcno = 0.10, X^ = 0.03, and Z = 0.02 (case Ne 2 in Table El. Each 
panel shows (a) V magn itude (open circles) taken from Chochol et al~lll993T) . (b) J, (c) H, and (d) K magnitudes, IR data (open circles) of which are taken from 
Woodward et al. 1 1997). One epoch is indicated by an arrow, i.e., when the optically thick winds stop for the 1.05 Mq WD model. Large open circles at the 
lower end of each free-free light curve indicate the epoch when the optically thick winds stop. A straight line of F\ oc r~ 3 is also added after the optically thick 
wind stops. 



iSala & Hernanz I (|2p05) calculated static sequences of hy- 
drogen shell-burning and compared the evolutional speed of 
post-wind phase for V1974 Cyg. They suggested that the 
WD mass is 0.9 M Q for 50% mixing of a solar composi- 
tion envelope with an O-Ne degenerate core (i.e., X = 0.35, 
Xcno = 0.25, and X Ne = 0.16), or 1 .0 M for 25% mixing (i.e., 
X = 0.53, Xcno = 0.13, and X^ e = 0.08). Their two cases re- 
semble our models of case Ne 1 and case Ne 2 in Table [2] 
respectively. Their WD masses are roughly consistent with 
our values. 

Several groups e stimated the WD mass of VI 974 Cyg. 
iRetter et al. I i ll 9971) obtained a mass range of Mwd = 0.75- 
1 .07 M Q based on the precessing disk model of superhump 
phenome non. A similar range of 0.75- 1.1 Mq is also ob- 
tained by Pare sce et al. I (119951) from various empirical rela- 



tions on novae. Both our 1.05 and 0.95 M© WD models are 
consistent with these constraints. 

5.5. Emergence of companion 

VI 974 Cyg is a binary system wit h an orbital period of 
P nrh = 0.0812585 d ays (1 .95 hr) (e.g.. iDe Young & Schmidt! 
1994: IRetter eTai~lll997l) . The companion mass is estimated 
to be M 2 = 0.15 M Q from equation d23l . For the 1.05 M 
WD model, the separation is a = 0.85 Rq and the effective 
radii of the WD Roche lobe and the secondary Roche lobe are 
R*\ = 0.44 R G and R% = 0.22 R Q , respectively. In the 0.95 M Q 
WD model, a = 0.82 R Q , R* = 0.44 R Q , and R* = 0.19 fl . 
The companion emerges from the WD envelope when the 
photosphere shrinks to ~ 0.8 Rq. This epoch is estimated 
to be day 95 and day 110 for the WD mass of 1 .05 Mq (case 
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FIG. 23. — Same as Fig |19l but for a different chemical composition of 
X = 0.35, X C no = 0.20, X Nc = 0. 10, and Z = 0.02 (case Ne 1 in Table|2|. The 
best-fit model is the 0.95 Mq WD among the three 0.9, 0.95, and 1.0 M Q 
WDs. 
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X = 0.35, Xcno = 0.20, X Ne = 0.10, and Z = 0.02 (case Ne 1 in Table|2). The 
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Ne 2) and 0.95 M Q (case Ne 1), respectively. 

ROSAT observations show that hard X-ray flux increases 
on day 70-100 and then decays on day 270-300. This hard 
component is sugge sted to have o ri ginated in the shock be- 
tween ejecta (Krau tter et"ai~ll 19961) . lHachisu & Kato I (120051) 
suggested a different idea that these hard X-rays originated 
from a shock between the optically thick wind and the com- 
panion. Before day ~ 70, the companion resides deep inside 
the WD photosphere and we do probably not detect hard X- 
rays. After the companion emerges, the shock front can be 
directly observed. Then the increase in the hard X-ray flux 
may correspond to the appearance of the binary component. 
The decrease in the hard X-ray flux may be caused by the 
decay of optically thick winds on day ~ 270-280. 

6. DISCUSSION 

6.1. Difference in free-free constant 

In a nova explosion theory, the evolution of a nova depends 
on four parameters, i.e., the WD mass, chemical composition 
of the envelope, mass accretion rate from the companion star, 
and thermal condition of the WD before the ignition. Our 



light-curve model basically follows a nova evolution after the 
nova envelope settles down in a steady state, i.e., after some 
time has passed from the optical peak. Therefore, the main pa- 
rameters that govern the evolution of novae are reduced to two 
from four, i.e., the WD mass and chemical composition. The 
other two parameters play an important role in the very early 
phase but do not affect the evolution solution in the steady- 
state phase. These two parameters are closely linked with the 
ignition mass and govern the free-free emission parameter, C\ 
in eauationfTot. because it is determined by the optically thin 
ejecta outside the photosphere. This is the reason why C\ val- 
ues are different among novae, and this is a subject of a new 
project. 

6.2. CNO abundance 

The CNO abundance is taken as a whole because the in- 
dividual ratios of C, N, and O do not affect the evolution of 
novae if the total amount of CNO is unchanged. It is because 
the energy generation of the CNO cycle depends on the to- 
tal amount of CNO but hardly on the individual ratios in the 
steady-state phase of novae. The CNO cycle changes only the 
relative ratio of each C, N, and O but keeps the total amount 
of C+N+O unchanged. Moreover, the Rosseland mean opac- 
ity is hardly changed if we choose another set of C, N, and 
O keeping the total amount of C+N+O constant. Therefore, 
the individual ratios of C, N, and O hardly affect the envelope 
evolution. 

We assume that, after the optical peak, the chemical com- 
position is constant throughout the envelope (everywhere) and 
throughout the nova evolution (everytime). We expect that C 
or O (or Ne) is dredged up from the WD interior in the very 
early phase of nova outbursts and then mix ed into the entire 
envelope by convection (e.g. Prialnik 1986). Our assumption 
means that the convection descends quickly after the optical 
peak and that the envelope becomes radiative in the decay 
phase of novae, in which the processed helium concentrates 
under the hydrogen-burning zone. 

From the computational point of view, we have some re- 
strictions in the chemical compositions of OPAL opacity. 
Therefore, we assume the total amount of CNO as a whole 
set in order to reduce our computer time and have selected 
several sets of chemical compositions as shown in table|2] 

7. CONCLUSIONS 

We propose a light curve model based on free-free emission 
and on the optically thick wind model, and have successfully 
applied it to three well observed novae, VI 500 Cyg, VI 668 
Cyg, and V1974 Cyg. Our main results are summarized as 
follows: 

(1) We have calculated light curves of novae in which free- 
free emission from optically thin ejecta dominates the contin- 
uum flux. The free-free luminosity is obtained from the den- 
sity structures outside the photosphere, which are calculated 
using the optically thick wind model. 

(2) The free-free emission light curves are homologous 
among various white dwarf (WD) masses, chemical compo- 
sitions, and wavelengths (in optical and infrared). Therefore, 
we are able to represent a wide range of nova light curves by 
a single template light curve on the basis of one-parameter 
family, i.e., a time scaling factor. 

(3) The template light curve declines as F oc t 5 in the 
middle part (from ~ 2 to ~ 6 mag below the optical max- 
imum) and then as F oc f" 3 5 in the later part (from ~ 6 to 
~ 10 mag), where t is the time from the outburst in units of 
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TABLE 11 
Physical properties of three novae 



subj ec t/ obj ect 


units 


V IDUU L.yg 


V looo Lyg 


v iy /4 cyg 


v i y / 4 uy g 


outburst year 


year 


1975 


1978 


1992 


i 


outburst day 


JD 


2,442,653.0 


2,443,759.0 


2,448,674.5 





nova speed class" 1 




very fast 


fast 


fast 




h d 


days 


.. 2.9 


12.2 


16 


, 


f 3 a 


days 


3.6 


24.3 


42 





early osc. 




no 


no 


no 





transition osc. 




no 


no 


no 





dust 




no 


thin dust 


no 




orbital period 


hours 


.. 3.35064 


3.32 


1.9488 




secondary mass^ 




0.29 


0.29 


0.15 




obs. WD mass 


Mq 


.. > 0.9 




0.75-1.1 


t 


E(B-V) 


.. 0.45 


0.40 


0.32 




obs. distance 


kpc 


.. 1.5 




1.8-1.9 




distance from UV fitting 


kpc 




3.6 


1.7 


1.6 


'break of y-band 


day 


.. 70 


86 






cal. WD mass 


Mq . 


.. 1.15 


0.95 


1.05 


0.95 


wind phase 


days 


.. 180 


280 


280 


260 


H-burning phase 


days 


.. 380 


720 


720 


590 


separation 


R e 


.. 1.28 


1.21 


0.85 


0.82 


companion's emergence 


days 


.. 50 


100 


95 


110 


chemical composition 




.. Ne2 


CO 3 


Ne2 


Ne 1 


hydrogen content (X) 




.. 0.55 


0.45 


0.55 


0.35 



a taken from lWamer] > 1995) for V1500 Cyg and V1974 Cyg but fromlMallama & Skillmanl H9791 for V1668 
Cyg 

k estimated from equation 1231 
see Tabled 



days. This break on the light curve is caused by the sharp de- 
crease in the wind mass loss rate. Since the time of the break 
is proportional to the time scaling factor, we use this time of 
break to specify the timescale of a nova light curve in the one- 
parameter family instead of the time scaling factor. 

(4) If we know, from observation, the time of break in 
the light curve, we can tell when the optically thick winds 
stop (supersoft X-ray turn-on time) and when hydrogen shell- 
burning ends (supersoft X-ray turnoff time). We can also tell 
the duration of the ultraviolet (UV) burst phase. These charac- 
teristic timescales are uniquely specified by the time of break. 

(5) The empirical observational formula between ?2 and tj, 
i.e., ? 3 = (1.68 ±0.08) r 2 + (1.9± 1.5) days, is derived from a 
slope of F oc t" 1 . 

(6) Our modeled light curves have been applied to three 
well-observed novae, VI 500 Cyg (Nova Cygni 1975), VI 668 
Cyg (Nova Cygni 1978), and V1974 Cyg (Nova Cygni 1992). 
Direct fittings with our free-free light curves indicate the WD 
mass of M WD « 1 . 15 M for V1500 Cyg and M WD « 0.95 M 
for VI 668 Cyg for the chemical compositions suggested. In 
V1974 Cyg, M WD w 1.05 M Q is obtained for the hydrogen 
content of X = 0.55 while M WD « 0.95 M Q is suggested for 
the different hydrogen content of X = 0.35. 

(7) Model light curves of the 1455 A band nicely follow 
the observations both for V1668 Cyg and V1974 Cyg. Fitting 
with the UV 1455 A light curves, we estimate the distances 



of V1668 Cyg to be d = 3.6 kpc for E(B-V) = 0.40 and the 
distance of V1974 Cyg to be d = 1 .7 kpc for E(B-V) = 0.32. 

(8) The supersoft X-ray flux was observed in V1974 Cyg, 
which emerged on day ~ 250 and declined on day ~ 600. 
This feature is explained consistently by our model. The pho- 
tospheric temperature rises high enough to emit supersoft X- 
rays on day ~ 250 - 280, which corresponds to the end of 
optically thick winds. The X-ray flux keeps a constant peak 
value for ~ 300 days followed by a quick decay on day ~ 600, 
which corresponds to the decay of hydrogen shell-burning. 

(9) Finally, we strongly recommend observations with the 
medium band Stromgren y-filter to detect the break of the light 
curve because the y-filter cuts the notorious emission lines in 
the nebular phase and reasonably follow the continuum flux 
of novae. 



We thank A. Cassatella for providing us with the machine 
readable UV 1455 A data of V1668 Cyg and V1974 Cyg, 
and the American Association of Variable Star Observers 
(AAVSO) for the visual data of VI 500 Cyg, VI 668 Cyg, and 
V1974 Cyg. We are also grateful to the referee, Marina Orio, 
for the useful comments that improved the manuscript. This 
research has been supported in part by the Grant-in-Aid for 
Scientific Research (1654021 1, 16540219) of the Japan Soci- 
ety for the Promotion of Science. 



REFERENCES 



Allen, C. W. 1981, Astrophysical Quantities (3rd ed.; London: Athlone), 103 
Ando, H., & Yamashita, Y. 1976, PASJ, 28, 171 

Andreae, J., & Drechsel, H. 1990, Physics of Classical Novae, eds. A. 

Cassatella and R. Viotti (Berin: Springer- Verlag), 204 
Andrea, J., Drechsel, H., & Starrfield, S. 1994, A&A, 291, 869 
Arkhipova, V. P., & Zaitseva, G. V. 1976, Soviet Astron. Lett., 2, 35 
Arkhipova, V. P., Burlak, M. A., & Esipov, V. F. 2000, Astronomy Letters, 

26, 372 



Augusto, A., & Diaz, M. P. 2003, AJ, 125, 3349 

Austin, S. J., Wagner, R. M., Starrfield, S., Shore, S. N., Sonneborn, G., & 

Bertram, R. 1996, AJ, 111, 869 
Balman, S., Krautter, J., & Ogelman, H. 1998, ApJ, 499, 395 
Becker, H. J., & Duerbeck, H. W. 1980, PASP, 92, 792 
Bode, M. F. et al. 2006, ApJ, in press i astro-ph/0604618 i 
Capaccioli, M., della Valle, M., D'Onofrio, M., Rosino, L. 1990, ApJ, 360, 

63 



Decline Law of Classical Novae 



23 



Cassatella, A., Altamore, A., & Gonzalez-Riestra, R. 2002, A&A, 384, 1023 
Cassatella, A., Benvenuti, P., Clavel, J., Heck, A., Penston, M., Macchetto, 

R, & Selveffi, P. L. 1979, A&A, 74, L18 
Chochol, D., Grygar, J., Pribulla, T., Komzik, R., Hric, L., & Elkin, V. 1997, 

A&A, 318, 908 

Chochol, D., Hric, L., Urban, Z., Komzik, R., Grygar, J., & Papousek, J. 1993, 

A&A, 277, 103 
Cohen, J. G. 1985, ApJ, 292, 90 
Collins, P. 1992, IAU Ore., 5454 

de Freitas Pacheco, J. A., dell'Aglio Dias da Costa, R., & Codina, S. J. 1989, 
ApJ, 347, 483 

De Young, J. A., & Schmidt, R. E. 1994, ApJ, 431, L47 

Duerbeck, H. W., Rindermann, R., & Seitter, W. C. 1980, A&A, 81, 157 

Duerbeck, H. W., & Wolf, B. 1977, A&AS, 29, 297 

Ennis, D., Becklin, E. E., Beckwith, S., Elias, J., Gatley, I., Matthews, K., 

Neugebauer, G., & Willner, S. P. 1977, ApJ, 214, 478 
Ferland, G. J. 1977, ApJ, 215, 873 

Ferland, G. J., Lambert, D. L., & Woodman, J. H. 1986, ApJS, 60, 375 
Ferland, G. J., & Shields, G. A. 1978, ApJ, 226, 172 

Gallagher, J. S., Hege, E. K., Kopriva, D. A., Butcher, H. R., & Williams, R. 

E. 1980a, ApJ, 237, 55 
Gallagher, J. S., Kaler, J. B., Olson, E. C, Hartkopf, W. I., & Hunter, D. A. 

1980b, PASP, 92, 46 
Gallagher, J. S., & Ney, E. P. 1976, ApJ, 204, L35 

Gehrz, R. D., Hackwell, J. A., Grasdalen, G. I., Ney, E. P., Neugebauer, G., 

& Sellgren, K. 1980, ApJ, 239, 570 
Hachisu, I., & Kato, M. 2000a, ApJ, 536, L93 
Hachisu, I., & Kato, M. 2000b, ApJ, 540, 447 
Hachisu, I., & Kato, M. 2001a, ApJ, 553, L161 
Hachisu, I., & Kato, M. 2001b, ApJ, 558, 323 
Hachisu, I., & Kato, M. 2003a, ApJ, 588, 1003 
Hachisu, I., & Kato, M. 2003b, ApJ, 590, 445 
Hachisu, I., & Kato, M. 2003c, ApJ, 598, 527 
Hachisu, I., & Kato, M. 2004, ApJ, 612, L57 
Hachisu, I., & Kato, M. 2005, ApJ, 631, 1094 
Hachisu, I., & Kato, M. 2006a, ApJ, 642, L53 

Hachisu, I., Kato, M., Kato, T., & Matsumoto, K. 2000, ApJ, 528, L97 

Hachisu, I., Kato, M., & Nomoto, K. 1996, ApJ, 470, L97 

Hachisu, I., Kato, M., & Nomoto, K. 1999a, ApJ, 522, 487 

Hachisu, I., Kato, M., Nomoto, K, & Umeda, H. 1999b, ApJ, 519, 314 

Hachisu, I., Kato, M., & Schaefer, B. E. 2003, ApJ, 584, 1008 

Hassall, B. J. M., et al. 1990, Physics of Classical Novae, eds. A. Cassatella 

and R. Viotti (Berin: Springer- Verlag), 202 
Hauschildt, P. H, Starrfield, S., Shore, S. N., Allard, E, Baron, E. 1995, ApJ, 

447, 829 

Hayward, T. L., Saizar, P., Gehrz, R. D., Benjamin, R. A., Mason, C. G., 
Houck, J. R., Miles, J. W., Gull, G. E., & Schoenwald, J. 1996, ApJ, 469, 
854 

Home, K, & Schneider, D. P. 1989, ApJ, 343, 888 

Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 

Iijima, T. 2006, A&A, 45 1 , 563 

Kaluzny, J. 1990, MNRAS, 245, 547 

Kato, M. 1983, PASJ, 35, 507 

Kato, M. 1994, A&A, 281, L49 

Kato, M. 1997, ApJS, 113, 121 

Kato, M. 1999, PASJ, 51, 525 

Kato, M., & Hachisu, I., 1994, ApJ, 437, 802 

Kato, M., & Hachisu, I., 2005, ApJ, 633, LI 17 

Kawara, K, Maihara, T., Noguchi, K., Oda, N., Sato, S., Oishi, M., & Iijima, 

T. 1976, PASJ, 28, no. 1, 1976, p. 163 
Klare, G., Wolf, B., & Krautter, J. 1980, A&A, 89, 282 
Krautter, J., Ogelman, H., Starrfield, S., Wichmann, R., & Pfeffermann, E. 

1996, ApJ, .456, 788 
Lance, C. M., McCall, M. L., & Uomoto, A. K. 1988, ApJS, 66, 151 
Lockwood, G. W., & Millis, R. L. 1976, PASP, 88, 235 



Lyke, J. E. et al. 2001, AJ, 122, 3305 

Lyke, J. E. et al. 2003, AJ, 126, 993 

Mariama, A. D., & Skillman, D. R. 1979, PASP, 91, 99 

Morisset, C, & Pequignot, D. 1996, A&A, 312, 135 

Morrison, W. 1978, IAU Circ, 3264 

Osborne, J. et al. 2006a, ATel, 764 

Osborne, J. et al. 2006b, ATel, 838 

Orio, M. 2004, Revista Mexicana de Astronorma y Astroffsica (Serie de 

Conferencias), Vol. 20. IAU Colloquium 194, pp. 182-186 
Paresce, E, Livio, M., Hack, W., & Korista, K. 1995, A&A, 299, 823 
Petitjean, P., Boisson, C, & Pequignot, D. 1990, A&A, 240, 433 
Pfau, W. 1976, Inf. Bui. Variable Stars, 1106 
Prialnik, D. 1986, ApJ, 310, 222 
Prialnik, D., & Kovetz, A. 1995, ApJ, 445, 789 
Retter, A., Leibowitz, E. M., & Ofek, E. O. 1997, MNRAS, 286, 745 
Rieke, G. H, & Lebofsky, M. J. 1985, ApJ, 288, 618 

Saizar, P., Starrfield, S., Ferland, G. J., Wagner, R. M., Truran, J. W., Kenyon, 
S. J., Sparks, W. M., Williams, R. E., & Stryker, L. L. 1991, ApJ, 367, 310 

Saizar, P., Starrfield, S., Ferland, G. J., Wagner, R. M., Truran, J. W., Kenyon, 
S. J., Sparks, W. M., Williams, R. E., & Stryker, L. L. 1992, ApJ, 398, 651 

Saizar, P., Pachoulakis, I., Shore, S. N., Starrfield, S., Williams, R. E., 
Rothschild, E., & Sonneborn, G. 1996, MNRAS, 279, 280 

Sala, G., & Hernanz, M. 2005, A&A, 439, 1057 

Schmidt, Th. 1957, Z. Astrophys., 41, 181 

Schwarz, G. J. 2002, ApJ, 577, 940 

Schwarz, G. J., Starrfield, S., Shore, S. N., Hauschildt, P. H. 1997, MNRAS, 
290, 75 

Scott, A. D., Duerbeck, H. W., Evans, A., Chen, A.-L., de Martino, D., 
Hjellming, R., Krautter, J., Laney, D., Parker, Q. A., Rawlings, J. M. C, 
& Van Winckel, H. 1995, A&A, 296, 439 

Seaton, M. J. 1979, MNRAS, 187, 73 

Semeniuk, I., Olech, A., & Nalezyty, M. 1995, Acta Astron., 
Shanley, L., Ogelman, H, Gallagher, J. S., Orio, M., & Krautter, J. 1995, ApJ, 
438, L95 

Shore, S. N. et al. 2003, AJ, 125, 1507 

Slavin, A. J., O'Brien, T. J., Dunlop, J. S. 1995, MNRAS, 276, 353 
Slovak, M. H. & Vogt, S. S. 1979, Nature, 277, 114 

Snijders, M. A. J., Batt, T. J., Roche, P. E, Seaton, M. J., Morton, D. C, 
Spoelstra, T. A. T., & Blades, J. C. 1987, MNRAS, 228, 329 

Stickland, D. J., Penn, C. J., Seaton, M. J., Snijders, M. A. J., & Storey, P. J. 
1981, MNRAS, 197, 107 

Tempesti, P. 1979, Astronomische Nachrichten, 300, 51 

Tomkin, J., Lambert, D. L., & Woodman, J. 1976, A&A, 48, 319 

Tylenda, R. 1978, Acta Astron., 28, 333 

Umeda, H, Nomoto, K, Yamaoka, H, & Wanajo, S. 1999, ApJ, 513, 861 
Vanlandingham, K. M., Starrfield, S., & Shore, S. N. 1997, MNRAS, 290, 87 
Vanlandingham, K. M., Starrfield, S., Shore, S. N., Sonneborn, G. 1999, 
MNRAS, 308, 577 

Vanlandingham, K. M., Schwarz, G.J., Shore, S.N., Starrfield, S., & Wagner, 

R.M. 2005, ApJ, 624, 914 
Wade, R. A., Ciardullo, R., Jacoby, G. H, & Sharp, N. A. 1991, AJ, 102, 

1738 

Warner, B. 1995, Cataclysmic variable stars, Cambridge, Cambridge 

University Press 
Williams, R. E., & Gallagher, J. S. 1979, ApJ, 228, 482 
Williams, R. E., Ney, E. P., Sparks, W. M., Starrfield, S. G., Wyckoff, S., & 

Truran, J. W. 1985, MNRAS, 212, 753 
Williams, R. E., Woolf, N. J., Hege, E. K, Moore, R. L., & Kopriva, D. A. 

1978, ApJ, 224, 171 
Woodward, C. E., Gehrz, R. D., Jones, T. J., Lawrence, G. F, & Skrutskie, 

M. F. 1997, ApJ, 477, 817 
Young, P. J., Corwin, H. G., Bryan, J., & de Vaucouleurs, G. 1976, ApJ, 209, 

882 



